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ABSTRACT 


The satellite to satellite tracking data unique processing considerations are 
presented. The strategy to be employed in the application of data processing 
routines to the ATS-F/GEOS-C satellite to satellite tracking (SST) experiment 
is described. When these considerations are included, studies, analysis and 
simulations predict an orbital position uncertainty of less than 15 meters and 
a velocity uncertainty of less than 1.5 millimeters per second. 
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INTRODUCTION 


The Satellite to Satellite Tracking experiment to be performed during the 
GEOS-C/ATS-F mission is primarily designed to produce data which can be 
utilized to derive coefficients for the harmonic expansion of the gravity field 
and in a related effort to aid in the GEOS-C program for Geoid shape deter- 
mination using the altimeter data/ Both of these goals or objectives require 
highly accurate orbits over long arcs? Since GEOS-C will be visible from 
ATS-F for over one half an orbit these long data spans are realistic. 

The geometry for this satellite to satellite tracking experiment is shown 
in Figure 1. 



Figure 1. GEOS-C/ ATS-F/GS Geometry 


The ATS-F satellite will be in synchronous orbit over 94° West longitude during 
the early portion of the GEOS-C mission. Later ATS-F will be drifted to 34° 
East longitude. It is planned to track GEOS from ATS-F from both positions. 

The nominal GEOS-C orbit will be;^ 

Height =' 843 km 
Eccentricity = 0.004 
Inclination = 115° 

The portion of the GEOS orbit visible from ATS-F will vary with time as the two 
satellites revolve in their respective orbital planes. 
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SPACECRAFT SYSTEMS 


The transponders in the two spacecrafts and the associated antennas establish 
the tracking capability of this experiment. 

ATS-F System 


The ATS-F spacecraft^ will employ a 9 meter (36 dB gain, 1.4° beamwidth) 
parabolic antenna programmed to follow the line of sight to the GEOS-C spacecraft. 
The communication transponder aboard the ATS-F will coherently translate the 
ground to ATS-F uplink frequency of approximately 6 GHz to 2 GHz for the 
ATS-F to GEOS forward link. On the return link the ATS-F transponder trans- 
lates the GEOS-C to ATS-F return link from 2.25 GHz to 4 GHz. The 4 GHz 
signal is the ATS-F to ground station link frequency. These frequency trans- 
lations are termed "pseudocoherent" in that the forward link translation is 
accomplished employing a phase locked oscillator that is coherent with the 
forward link signal at the time of translation, and the return link translation is 
accomplished employing the same phase locked oscillator however the oscillator 
is not coherent with the return link signal. The ATS-F translation oscillator 
is always phase locked to the ground source during satellite to satellite tracking 
operations. 

GEOS-C System 

The GEOS-C tracking transponder is a phase locked frequency multiplication 
device with a 240/221 transmit to receive ratio. The GEOS-C is fitted with four 
switchable antenna on the slant panels for receiving and transmitting to ATS-F. 

A fifth S-band antenna is moimted on the earth viewing side of GEQS-C. When 
commanded to this antenna the transponder may be interrogated directly from 
the ground. 

The tracking data recorded at the Ground Station is in the form of time 
delays and doppler frequency count times corresponding to the roimd trip range 
and the roimd trip range rate (Doppler). 

PREPROCESSING CONSTRAINTS 

Each raw satellite to satellite tracking data record consists of two basic 
time measurements. These are usually referred to as range and range rate. 

The actual "range” measurement is a measure of the two way signal propagation 
time including systematic or instrument delays of a specific event, such as the 
zero crossing of a range tone. The actual "range rate" measurement is a 
measure of the time required to accumulate a known number of doppler plus 
bias frequency cycles (N-count). Since the Application Technology Satellite (ATS) 
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ground station employs the N-count system the processing is designed for 
N-count doppler. 

Actual tracking geometry (Figure 1) shows that the observed "range” (time 
delay) is a function of the positions of the two satellites and the ground station. 

It is not possible to separate any individual satellite to satellite or satellite to 
ground propagation time. The processing approach is to treat the observation 
"as is" and let the orbit determination program synthesize a propagation time 
based upon a priori or ephemeris derived information. The same is true of the 
doppler count time and it is also treated "as is" in the processing. The data 
processor is configured to allow independent time tagging of the two observables. 

PREPROCESSING FUNCTIONS 

The preprocessing performs the following functions; 

(a) Read and correct data format 

(b) Provide proper time tag 

(c) Apply known corrections 

(d) Edit, smooth, and compact data 

(e) Format output data. 

The preprocessor approach to each of those is explained in the following sections. 

A. Input 

The fimction accomplished here is to read in the raw data in the ATSR data 
format and convert it to the binary Unified Data Format. 

B. Time Tagging 

The time tag associated with the raw data is the ground clock time for the 
transmission of a positive going zero crossing of a range tone. This same time 
is associated with the doppler since the doppler count period starts also at 
this zero crossing time. The time tagging algorithm associated with Satellite 
to Satellite tracking must take into accovint all the motions of the two spacecrafts 
and the ground station. This time tagging algorithm is reflected in the observational 
equations. The observational equation is simply an expression relating the 
observed measurement to a computed measurement. The computed measurement 
is a function of the two satellite positions at the time the signal arrives at each 
satellite, the ground station position (coordinates) and other system parameters. 

The geometry is illustrated in Figure 2 where 

Xj (T) = ATS position vector 
x^ (T) = GEOS-C position vector 

= Ground Station position vector. 
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Figure 2. Position Vectors 


The range tone is transmitted from the ground station at time (t) and reaches 
ATS-F at (t + Ti). Thus Tj represent the propagation time from the ground 
station to ATS-F. After a transponder delay of TDj , the signal is transmitted 
to GEOS-C arriving at t + + TD^ + Tj . Where Tj is the propagation time 

from ATS-F to GEOS-C. If the GEOS-C transponder delay is TDj , the signal 
arrives back at ATS-F at t + Tj + TD^ + Tj + TDj + T3 . Since the ATS-F 
transponder may not be exactly reciprocal the time delay TD3 is assigned for 
the return link delay. The signal then arrives back at the ground station at 
t + Tj + TDj + Tj + TD^ T3 + TD3 + . One more time delay (Tj^ ) is defined 

as the propagation delay caused by transmission path anomalies such as the 
troposphere. These time delays may be interpreted in terms of the various 
position vectors of Figure 2 as: 
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( 1 ) 


Tj = [XjCt + Tj) - Xj/c + Tjj 

Tj = [X^ (t + + TD^ + Tj) - X^(t + Tj + TDi)} /c 

T3 = [Xi (t + Tj + TDj + T2 + TD2 + T3) - XjCt + Tj + TD3 + + TD2)] /c 

= [X3 (t + Tj + TD^ + + TDj + T3 + TD3) - Xj /c + 

The solution of this system of nonlinear equation may be obtained by means of the 

Newton-Ralphson iterative method for the solution of nonlinear equations. 

Models for the propagation path anomalies must also be included to evaluate Tr . 

Data Corrections 

The orbit dependent variation in the transponder delays TD^ , TDj , TD3 are 
functions of the deviation of the signal frequency from its nominal value. Such 
variations in carrier frequency are due to the relative motion between trans- 
mitter and receiver (Doppler shift). These variations can be expressed as 
function of the time derivative of the propagation time delay. The accepted 
expression for average doppler is: 


f 


<1 


2ff. 


c 


( 2 ) 


Since the Doppler shift is a linear function of the transmission frequency, 
may be defined to be the sum of two frequencies; an fj which propagates 
from the ground to ATS-F to GEOS-C and return and anfj which propagates 
from the ground to ATS-F and return. The values assigned to and £2 and the 
proper frequency formulations are explained in the "Doppler Processing" sec- 
tion of this report. The Doppler expressed in terms of fj , fj and the propaga- 
tion paths of Figure 2 (r, and r^) is 

fd = 

If we now define the time derivative of the propagation delay as T^. 

fd " - + Tj) + ^2 (ti)] ( 4 ) 


( 3 ) 

f. 

= Vc then 
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Since the ATS-F is in geos 5 mchronous orbit most of the observed Doppler is 
from the ATS-F/GEOS-C link. For evaluation of the variations of the various 
transponder delays due to the Doppler, the Ground to ATS-F Doppler will be assumed 
to be zero. However since, as stated before, it is impossible to separate the 
various contributions equation 4 for one way Doppler is approximated as; 

f, = fi[Ti4T2] (5) 

As previously stated the ATSR system employs ’'N-count” Doppler. The 
system determines the time, (T), required to accumulate a fixed number of 
Doppler plus bias frequency cycles (N). That is 


T = 


N 


fd 


v/hich when combined with equation 5 gives: 

N 

which gives an approximate doppler count of 


( 6 ) 


and 


T = 


N 


f 2R , 


(7) 


K = 2^(N-Tfj,) (8) 

This expression for R is sufficiently accurate for the determination of 
delay variation in the ATS-F transponder due to Doppler frequency offset. It 
is to be noted that the dynamic loop lag of the GEOS-C transponder must also be 
in the modeling TD^ . The transponder delay variations can be given by the 
expressions: 
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( 9 ) 


4 



4 ^ 


k = 0 
4 


TD. 




k -0 


where d,^ and e^ are a set of polynomial coefficients associated with the ATS-F 
and the GEOS-C transponder and represents the GEOS^C phase lock loop 
acceleration gain constant. The time derivatives {T. ) may be evaluated from 
the following expressions : 

Tj = [Xj(t + T,) - XJ • Xj (t + Tj):/c2Tj 

= [X^(t^ + T,) - Xi(t^)3 • [X^Ct^ + Tj) -X,(t^)]./c2T2 

T3 = [Xj(t^ +T2) -X^ct^)] ■ [Xi(t^ +X,)-X2(t,)]/c2T3 

T4 = tXi(t,)-Xj -X^Ct^Vc^T, 

■ftiiere 

tq = t + Tj + TDj 
= ,t^ +T3 + TD3 

Thus the range observation, expressed in time delays is; 

T = 13 + TDj + Tj + TDj + T3 + TD3 + T4 (11) 

The lag error is a function of the range acceleration and the loop gain constant. 
The loop gain constant is of course a loop design parameter and is therefore 
readily available. However the range acceleration must somehow be developed 
either from the data or the orbit. The proposed method of obtaining the range 
acceleration is outlined in the following section. 
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The data corrections derived from the station and instrument calibrations 
must also be applied to the observational data. These corrections include: the 
station survey corrections in East, North and Vertical components, instrument 
bias, timing; bias and the previously .mentioned refraction bias. These corrections 
are modeled in the NAP-3 and are applied to the appropriate measurements by 
the program. The value of t is transmitted in the ATSR data format'^ as Words 
2 and 5 in nanoseconds. 

SMOOTHING AND EDITING 

The routine employed for data smoothing is a portion of the Navigation 
Analysis Program (NAP).’ The data are fitted to a polsmomial function of 
time. The number of data points to be included in the fitting as well as the order 
of the polynomial (up to 6th order in t) may be selected as user inputs. The user 
may also select the criteria for the rejection of statistically inconsistant data 
points . 

Once this polynomial fitting is accomplished the range acceleration may 
be derived directly since the acceleration is the slope of the fitted curve. As 
indicated in the previous section of this paper the range acceleration is required 
to evaluate the dynamic loop lag in the GEOS-C transponder. In utilizing this 
capability the user must select: 

(1) the order of the pol 3 momial 

(2) rejection criteria for statistically inconsistant data points 

(3) the number of data points to be smoothed 

With this acceleration model the transponder delay variation for the GEOS-C 
can be evaluated. 

DOPPLER PROCESSING 

The Doppler or time delay rate measurement analysis is similar but more 
involved. The propagation time delay is not effected by carrier frequency 
translations and multiplications. However, since the Doppler is proportional 
to the carrier frequency the various carrier frequency adjustments must be 
considered. The frequency of the carrier may be traced with the aid of 
Figure 3. The signalsappearing at the various points in the diagram are listed 
in Table 1. The constants involved are given in Table 2. The signal at x (Table 
1) is the input to the phase locked ATSR receiver for doppler extraction. The 
extracted doppler is then, 

= (^)(^^) + K3) + b(K2 - Kj) + c(Kj + K3)] - 1 ] 

( 12 ) 
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TABLE 1 


A (K^ + K3) f, 

B (Ki + K 3 ) f, 

C ai(K2 + K3)f^ 

D (K, + K3) 


E 

F 


ba,(K, + K3) 
^'' 3 ) 

bajCK^+Kg) 


aj(K^ + K3) ^ ^ ^ j ^ ^ 

^(Ki +K 3 ) 7 


fr = [2(^2 + K3) + b(K2 - Kj)] 

3 D 

fr = + Ka)^ fr 


K 


a + b 
ca,(K^ + Kj) 


a + b 


+ — ^- [a(K2 + K3) + b(K2 - Kj)] 7 = 

3 “|- D 


[a(K2 + K3) + bCKj - K,) + c(K, + K3)] f^= f, 

a /Si 

L / 3 j 7 =J-Jl [a(K2 + K3) + b(K2 - K^) + c(K, + I^)] 
a + b 


N 

^ m 


R " f, 

^ m ^ 






m 


W + d - e aj (Ki + K3) f^Oj 


m 


“ 2 AA t (7i) “■ A - A 
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TABLE 2 
f = Sx 10^ Hz 

r 

K, = 16.43 
K 2 = 15.2525 
K 3 = 1213.57 
= 775.4 
a = 6 

b ^ 240 
c = 77 
d=152 
e = 84 
m = 221 
n = 240 


- 1 tj 


<^2 = 


A - 1 - T3 


reference 5 




1 + T. 


V J 


11 



In the implementation of the system a bias frequency is added to this Doppler 
frequency to prevent the count from going through zero. If the count goes 
throxigh zero the direction of the range rate is lost. Inserting the constants 
of Table 2 in equation 12; 

= 22A1 - 1] + 1700[aja2 - 1] MHz ( 6 ) 

In the system (ATSR) data output the Doppler data appears as the time 
required to accumulate a known number of bias (fj^) plus Doppler frequency 
cycles. The output is 


T = 


N 


fb + fd 


where 

f,^ = bias frequency (500 KHz) 
f^ = Doppler frequency (equation 6 ) 

N = cycles counted (reference 4) 

In the implementation of this experiment the values of the constants Ki , Kj , K 3 
and K. cannot be realized exactly. Therefore a delta in the bias frequency must 
be introduced into the data system. The bias frequency is, 


where 


and 


f^ = correct doppler bias including 
frequency offsets 

f^ = nominal bias (500 KHz) 

A f = bias offset caused by transmitted 
and local oscillator offsets 



(14) 


12 



where 


A = frequency offset in the ranging carrier 
from ground to ATS-F 

A f^^g = frequency offset in the ground to ATS-F 
beacon 

A = frequency offset in the first down 
converter in the ATSR receiver 

These frequency offsets differ with different operating frequencies and also with 
the ground stations. As of this time there are no plans to include these as part 
of the data message. The computer program is designed to accept the A f^ as an 
input variable. 

For the selected frequencies in the ATS-F/GEOS-C experiment this Af^^ = 
+16.3102 and the zero Doppler count, T - .51190330 seconds. 

CONCLUSION 

The test of any system is ”How well can the desired mission be accomplished?" 
To answer this question, the orbit determination capability of the system was 
investigated employing the NAP-3 (ref. 7). This program was developed to 
accept and process various types of tracking data. As such is includes the 
algorithms required to process satellite to satellite range sum and range 
sum rate data. The NAP-3 actually considers the four paths separately and 
through an iterative solution resolves the final orbit for each of the two satellites. 
This capability has been exercised and evaluated using simulated data.® 

Studies® involved in these simulations have indicated the necessity of solving 
not only for the two satellites but also the more dominant terms of the harmonic 
expansion of the geopotential. The question then becomes, ^'What is the smallest 
number of geopotential harmonic expansion coefficients to be considered to 
meet a position uncertainty criteria?" A position uncertainty of 25 meters 
was chosen for our analysis. A covariance analysis (ref. 10) was employed 
to answer this question. The object of this analysis is to satisfy the constraint 
of 25 meters total position imcertainty by adding the smallest possible number 
of coefficients to the adjusted parameter set which will accomplish the task. 

The procedure employed, as described in reference 11 is recursive where 
the (N + 1)®’^ coefficient is adjusted by examining the results of adjusting N 
coefficients along with the GEOS-C and ATS-F states. 

When gravity harmonics are included in the solution, position and velocity 
uncertainties of Figures 4 and 5 are acheived. 
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